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■ ABSTRACT 

czj ! 

C3 . Context. The runaway star HD 34078, initially selected to investigate small scale structure in a foreground diffuse cloud has been shown to be 
surrounded by highly excited H2, the origin of which is unclear. 
I Aims. We first search for an association between the foreground cloud and HD 34078. Second, we extend previous investigations of temporal 
• absorption line variations (CH, CH + , H2) 
j-H ' Methods. We have mapped the 12 CO(2 — 1) emission at 12" resolution around HD 34078's position, using the 30 m IRAM antenna. The 
, follow-up of CH and CH + absorption lines has been extended over 5 more years. In parallel, CH absorption towards the reddened star £ Per 
' have been monitored to check the homogeneity of our measurements. Three more FUSE spectra have been obtained to search for NQI2) vari- 
es ations. 

^v^l . Results. CO observations show a pronounced maximum near HD 34078's position, clearly indicating that the star and diffuse cloud are as- 
f"*) ' sociated. The optical spectra confirm the reality of strong, rapid and correlated CH and CH + fluctuations. On the other hand, 7V(H2, J = 0) 
^\ has varied by less than 5% over 4 years, indicating the absence of marked density structure at scales below 100 AU. We also discard iV(CH) 

■ variations towards £ Per at scales less than 20 AU. 

^ Conclusions. Observational constraints from this work and from 24 fj,m dust emission appear to be consistent with H2 excitation but incon- 
. , sistent with steady-state bow shock models and rather suggest that the shell of compressed gas surrounding HD 34078 or lying at the boundary 
1 of a small foreground clump, is seen at an early stage of the interaction. The CH and CH + time variations as well as their large abundances 
^ ' are likely due to chemical structure in the shocked gas layer located at the stellar wind/ambient cloud interface. Finally, the lack of variations 
for both NQI2, J = 0) towards HD 34078 and iV(CH) towards £ Per suggests that quiescent molecular gas is not subject to pronounced 
small-scale structure. 

Key words. ISM: molecules - stars: individual (HD 34078) - ISM: structure 



1. Introduction 

During the past decade strong evidence has accumulated indi- 
cating that the spatial distribution of species like Na I and Ca II 
within neutral interstellar (IS) gas displays significant structure 



* Based on observations made mainly at IRAM, Observatoire de 
Haute Provence (France), McDonald Observatory (USA) and with 
FUSE 



at AU scales dCrawfordll2003l: lLauroeschll2007l; |WeltyJ|2007; 
Weltv et al. 2008I) . Hiitself. within the cold neutral medium 



at least, shows such structure ( Frail et all 1 1994 iHeilesI 1 19971 
iDeshpandel 120071: IWeisberiT & Stani mirovid 120071) . In diffuse 
molecular gas, similar conclusions have been reached for trac- 
ers like H 2 CO, HCO+, and OH dMoore & Marscherl 1 1995 



Liszt & L ucas 2000). Whether or not these spatial variations 
correspond to true density structure [i.e. to local fluctuations 
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of n(H2)] is obviously of key importance for the modelling of 
physical and chemical processes within molecular gas. 

To investigate this question, a time variation study of H2, 
CH, and CH + interstellar absorption lines towards the 09.5V 
runaway star AE A ur, HD 34078 has b een undertaken by 
Rollinde et ail d2003l , hereafter R03) and|BoisseetalJ d2005j 
, hereafter B05). This bright star is significantly reddened [E(B- 
V) = 0.53] and its optical spect rum displays strong a bsorption 
lines from e.g. CH, CH+, CN dFederman et al.lll994l) . typical 
of diffuse molecular clouds. HD 34078 has a large proper mo- 
tion of /j, = 43 mas yr _1 , corresponding to a transverse velocity 
of 103 kms -1 or 22 AU yr _1 for a distance D — 530 pc 
(this value will be used in the following for consistency with 
B05; it is compatible with the trigonometric parallax estimate 



of 446^ 22 " pc). The line of sight is thus drifting rapidly through 
the foreground gas; successive column density measurements 
then provide a "cut" through the cloud revealing its spatial den- 
sity structure over scales which typically range from 1 to 100 
AU for time separations ranging from a few weeks to a few 
years. 

The five first FUSE spectra discussed in B05 showed that 
highly excited H2 gas is present along the line of sight, together 
with more standard quiescent gas at T ~ 80 K. The presence 
of significant amounts of H 2 with an excitation energy higher 
that 2500 K (corresponding to the v = 0, J = 5 state) is very 
rare (another rema rkable case is that of HD 37903 studied by 
Meyer et al. 200lb . This is certainly related to the fact that by 
chance, along its long path from th e Orion nebula where i t 
was ejected a few millions years ago (IHoogerwerf et al.ll200ll) . 
HD 34078 has recently encountered a dense interstellar cloud 
with which it is currently interact ing, leading to the present-day 
apparence of the IC 405 nebula dHerbiglll958l) . 

In previous modelling work of the properties of the gas 
along the line of sight (B05), it was found that two components 
are required to account for the observed absorption lines: 

- highly excited H2 located in a bow shock, where the stellar 
wind impacts onto the ambient medium. Illumination of this 
gas by the strong UV field of the star satisfactorily explains the 
observed H2 excitation diagram. The presence of such a bow 
shock around HD 34078 wa s suspected on the basis of IRAS 
data (Ivan Buren et al.ll 1995b and has been confir med recently 



by higher spatial resolution Spitzer observations ( France et al 
2007), 



- a foreground quiescent cloud (supposedly unrelated to 
HD 34078/IC 405) giving rise to absorption from cold H2 and 
other species typical of translucent material (CH, CO, CN). 

However, in this scenario, the close similarity of the ra- 
dial velocity of both components appears as a pure coinci- 
dence. Further, the CH/H 2 ratio towards HD 34078 displays an 
anomalously high value, about three to four times larger than 
the one commonly found for other lines of sight. This led us 
to suspect that HD 34078 might be more closely related to the 
cold cloud than assumed in the above model, even if at first 
sight, this appears difficult to reconcile with the low H2 tem- 
perature (T = 77 K) derived from the population of the J = 
and J = 1 levels and the presence of a variety of molecules that 
HD 34078 might easily photodissociate. 



The potential relation between HD 34078 and the cloud 
probed has important implications regarding the investigation 
of small-scale structure. Indeed, in the case of a real associ- 
ation, any mechanical or radiative interaction might signifi- 
cantly affect the initial structure. Further, the gas flow in the 
bow shock may be subject to instabilities which could lead to 
time variations of a different nature, not necessarily associated 
with spatial structure. 

Thus, one objective of the present study is to clarify the re- 
lation between HD 34078 and the cold cloud. To this aim, we 
have undertaken observations of CO emission in the field sur- 
rounding the O star. The second goal of our work is to extend 
the search for absorption line variations in ground-based spec- 
tra (CH and CH + mainly) and in FUSE spectra (H2) by adding 
recent data obtained after the initial studies by R03 and B05 
were completed. 

This paper is organised as follows. We first present high 
spatial resolution 12 CO(2 - 1) and 12 CO(1-0) observations 
performed at the 30 m IRAM telescope and their implications 
concerning the location and properties of the quiescent com- 
ponent (Sect. 2). In Sect. 3, we analyse a new series of optical 
spectra which allow us to study in detail the variations of CH 
and CH+ column densities [hereafter N(CH) and iV(CH+)] 
and velocity profiles. Time variations of H2 column densities 
[hereafter TVfHa)] are then discussed in the light of the ob- 
served CH and CH + variations (Sect. 4). In Sect. 5, we sum- 
marize the main observational results for £ Per and HD 34078, 
in particular for readers not interested in details concerning ob- 
servations. Then, we discuss the implications of our observa- 
tions in terms of processes related to the interaction between 
HD 34078 and the surrounding gas and small-scale structure 
in foreground unperturbed gas (Sect. 6). Finally, we summa- 
rize our conclusions and present some prospects concerning 
possible observational signatures of the future evolution of 
HD 34078's close environment. 

2. CO emission towards HD 34078 

12 CO emission can be used to investigate a possible connec- 
tion between HD 34078 and the molecules seen in absorp- 
tion. Indeed, in the two-component model presented in B05, 
CO molecules (as well as CH, CH + etc.) lie far in front of 
HD 34078 and if this picture were correct, the morphology of 
the CO emission should not correlate in any manner with the 
star position. Otherwise, any relation between the CO emission 
map and HD 34078's position would be a clear indication that 
the star is closely related to the foreground cloud. 

2. 1. IRAM-30m observations 

Emission of the rotational transitions of 12 CO was observed at 
the IRAM-30m telescope during three consecutive nights on 
February 11 12 and 13, 2004. U sing the HERA 3x3 multi- 
beam array dSchuster et al J 120041) . we mapped the 12 CO(2-l) 
emission in a 66 x 66 arcsec 2 region centred on HD 34078's 
position. The single side-band receiver temperature was in 
the range 120 to 180 K. Observations were done under good 
weather conditions with 2 mm water vapor and a zenith sky 
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Fig. 1. ^ 2 CO(2-l) spectra o verlaid on the MIPS 24 fim emis- 
sion of France et al.l (12007b (color scale; note that the IR map 
is saturated near the maximum at the (3,12) position). The CO 
LSR velocity scale is in the range 0-13 km s _1 and the antenna 
temperature scale ranges from —0.1 to 0.6 K. The observed po- 
sitions correspond to the center of each panel; the dotted line 
indicates a constant velocity Vlsr = 5-9 kms -1 , associated 
with the CH and CH+ lines (see below). HD 34078 is located 
at offsets (0,0); the direction of its proper motion is indicated 
by the arrow. 



opacity at 230 GHz T230 — 0.15, resulting in system tempera- 
tures in the range 250-350 K. Chopper-wheel calibration was 
done every 10-15 minutes. Pointing was checked frequently, 
ensuring an accuracy of 2". We used the VESPA autocorrela- 
tor as a spectrometer covering 160 km s^ 1 with a resolution 
power of 7 x 10 5 or Sv = 0.4 km s _1 (8v = 0.3 MHz; see 
Table Q] for CO observation parameter values). The map was 
done in raster mode to allow deep integration of about 40 min. 
to ensure a rms of 15 mK (antenna temperature scale) in each 
velocity channel. It consists in a regular grid of 12 x 12 po- 
sitions observed with a sampling of 6", and the data are thus 
only slightly undersampled (HPBW = 11.7" at 230 GHz). 
A 5-point cross centred at offsets (—3, 3) with 12" steps was 
observed simultaneously with the single-pixel receivers facil- 
ity, in the 12 CO(1-0) and (2-1) transitions (HPBW = 22" at 
1 15 GHz) to derive the excitation conditions of the gas. 

2.2. Results and analysis of the CO data 

The resulting map of the 12 CO(2-l) emission is shown in 
Fig. [T| in which the nearly fully-sampled spectra were pro- 
jected on a 6 " x 6 " grid centred on the (0,0) of fset. Thermal 
dust em ission has also been observed recently by iFrance et al 
d2007h and their MIPS 24 /im is overlaid by the CO spectra. 




Fig. 2. The integrated 12 CO(2-l) e mission (contours) o verlaid 
on the MIPS 24 /im emission of France et al.l d2007l) (color 



scale). The CO profiles have been integrated in the interval 2.5 
- 10 km s _1 ; the first contour level and spacing between suc- 
cessive contours is 0.189 K.km s _1 . 



We also present a map of the integrated 12 CO (2 — 1) inten- 
sity superimposed onto the MIPS 24 /im emission (Fig. |2]i as 
well as channel maps (Fig.0. The 12 CO(2 — 1) integrated in- 
tensity is clearly stronger around the star position (0,0). While 
the NE and S W region are nearly devoid of emission, the region 
of enhanced 12 CO(2 — 1) emission is resolved and displays 
an elongated morphology (in the NE-SW direction) and extent 
which are both relatively similar to those of the brightest part of 
the IR arc. The 12 CO(2 — 1) peak is apparently offset south- 
ward by about 9 " with respect to the 24 /im maximum but, 
given the accuracy of absolute positions in the IR map (~ 6 "; 
K. France, private communication) and 12 CO(2 — 1) emission 
(~ 2 "), it is not clear whether this offset is really significant. 

The profiles are double peaked close to the star (in an ap- 
proximate circle of diameter 20-30 "; note the blue and red 
components appearing on each side of the vertical dotted line 
at 5.9 km s _1 plotted in Fig . |T} . The spectra closest to the star 
display a kind of mirror symmetry about an axis coincident 
with HD 34078's path, with the blue emission line stronger 
to the E and the red one stronger to the W. This is apparent 
in the channel maps (Fig. [3), where the red component peaks 
to the W while the blue one peaks to the NE. The symmetry 
quickly disappears as the profiles become single peaked further 
away. Averaging spectra in the NW or SE areas clearly shows 
that weak wings are present there over the entire velocity range 
(VLsr= 3-10 km s _1 ) covered by the more intense emission 
seen in the central part. The systemic velocity of the ambient 
molecular gas is Vlsr — 6.5 km s _1 . 

We now consider the additional 12 CO(1-0) spectra (Fig.@]i 
to investigate the excitation of the CO gas at the (-3,3) posi- 
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Fig. 3. Channel maps of the 12 CO(2-l) emission averaged over 
an interval of 0.41 km s _1 (the value in the upper left corner 
is the center of the velocity interval considered). Each panel 
corresponds to the region shown in Fig. Q] The first contour 
level and spacing between successive contours is 3.5 mK. 

tion (at other positions, the 12 CO(1-0) spectra are of lower S/N 
due to a smaller integration time and are less appropriate to 
measure the excitation ratio). For this purpose, 12 CO(1-0) and 
12 CO(2 — 1) spectra are brought in the main beam temperature 
scale; the forward and beam efficiencies are given in Table [T] 
We next synthetize the 12 CO(2 — 1) emission over the (1-0) 
beam using the spectra obtained at adjacent positions [(-9,3), (- 
3,9), (3,3), (-3,-3)] with the appropriate weighting. The result- 
ing (2-l)/(l-0) integrated intensity ratio appears to be around 
1 .5, a large value compa red to that (~ 0.7) commonly observed 



Velocity (km/ 



Fig. 4. Emission spectra of 12 CO(1-0) (grey scale histogram) 
and 12 CO(2 - 1) (thick line) centred at the (-3, 3) off- 
set. 12 CO(2 — 1) spectra from adjacent positions were con- 
volved to synthetize the 12 CO(2 — 1) emission in the 22" 
(1-0) Gaussian beam (see text). Main beam temperatures are 
adopted. 



towards diffuse clouds (Falgarone et al. 1998; Liszt & Lucas 



1998; Petv et alj|2008h . Moreover, the (2-l)/(l-0) ratio is sig- 
nificantly different for the two main components around Vlsr 
~ 5 and 7 km s _1 for which we get values of about 1.8 and 1.3, 
respectively. 

The above results on the (-3,3) position can be used 
to set constraints on the physical conditions prevailing in 
the gas. With the help of a large velocity gradient model 
(Hil y -Brant & Falgarondl2007D . we find that acceptable solu- 
tions can be obtained for T > 12 K and that above T = 20 K, 
the gas density is well constrained and must lie in the range n = 
10 3 — 10 4 cm -3 (n is the ambient H number density). All so- 
lutions obtained correspond to optically thin emission. For in- 
stance, at T = 80 K, a value close to the temperature estimated 
for the dust by lFrance et al. ( 2007 ). we get n = 2 x 10 3 cm~ 3 
and iV(CO) = 3.7 x 10 14 cm" 2 at Vlsr = 5 km s" 1 and n = 
1.8 x 10 3 cm- 3 and N(CO) = 6.0 x 10 14 cm~ 2 for the Vlsr 
= 7 km s _1 component. Note that the above CO column densi- 
ties are comparable to the value iV(CO) = 5.7 x 10 1 4 cm 



inferred either from IUE data by Mc Lachlan & Nand y ( 1984) 



or from FUSE by Shef fer et al.l d2008t) for the gas in front of 
HD 34078. 

We end this section by concluding that our CO observations 
allow us to unambiguously answer the question which motived 
them: HD 34078 is indeed closely associated with molecular 
gas located in its immediate vicinity (which was previously 
assigned to a foreground quiescent cloud by B05). Moreover, 
the anomalous CO excitation observed, the large inferred gas 
density and peculiar velocity field very likely result from the 
interaction between the stellar wind/radiation and the ambient 
molecular material (in Sect. 6, we discuss in more detail the 
implications of the CO observations). 



3. Variation of CH and CH absorption lines 

3.1. Description of optical observations 

In the visible, we add 26 spectra to the data considered in 
R03. Twelve spectra were taken at OHP, eight at McDonald 
Observatory (hereafter McD) while three spectra were obtained 
at the Boyunsan Optical Astronomy Observatory (BOAO, 
South Korea), two at the Terskol Observatory and one at Calar 
Alto. Altogether, these observations probe the evolution of CH 
and CH + abundances between 1989 and 2008, with good sam- 
pling since 2000. In particular, our recent data cover the period 
during which the eight FUSE spectra were obtained well. The 
date of each observation and spectral resolution are given in 
Table |2] 
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To check in a direct way the consistency of measurements 
performed at different telescopes, we started in 2003 paral- 
lel observations of the bright star ( Per. The latter has been 
observed in nearly all OHP and McD runs, in addition to 
HD 34078 (cf Table©. This nearby reddened star [d = 400 pc, 
E(B-V) = 0.33] has a small proper motion of 10.2 mas yr _1 , 
corresponding to a transverse velocity of 4.1 AU yr _1 , much 
smaller than the value for HD 34078 (22 AU yr" 1 ). Thus, we 
expect in principle, much less variation due to structure in the 
foreground IS gas towards £ Per and absorption lines seen in 
the spectra of that star should be a good indication of the instru- 
mental stability and homogeneity of our measurem ents. Data 
from the literature dAllenl 1 1 994c ICrane et al] 1 19951) do show 
that W\ values for CH, CH+, CN, Cal, and Call are constant 
within errors for £ Per. 

OHP observations were done in servic e mode using the 
ELODIE spectrograph dBaranne et al . 1996), as in R03, except 
for the four latest runs which were performed with SOPHIE, 
the new spec trograph that now supersedes ELO DIE at the 1 .93 
m telescope dBouchy & The Sophie Teaml2 006). SOPHIE pro- 
vides an improved spectral resolution (R = 75 000) and sen- 
sitivity, as well as an extended wavelength range (including 
CH A3886, CH A3890, and CH+ A3957; as the blue CN lines 
are close to the blue edge of the spectra the S/N is too low for 
these features to be usable). Since SOPHIE has been optimized 
for the detection of extrasolar planets by radial velocity mea- 
surements, it provides an accurate wavelength scale (better than 
0.01 km s _1 ; this scale is relative to the barycentre of the solar 
system). Each observation consisted in 4 to 8 individual expo- 
sures totalling about 1 hr for HD 34078 and 10™ for C Per; 
these two targets were generally observed consecutively or, oc- 
casionally, during two sucessive nights. Spectra were extracted 
using the pipeline data reduction software. The latter were 



designed specifically for these spectrogr aphs dBaranne et al 



1996; Bo uchy & The Sophie Team 2006) and include all re- 
quired steps (in particular bias and flat field corrections, using 
appropriate exposures). For observations of bright stars such as 
£ Per and HD 34078, this procedure has been shown to work 
efficiently, which we checked whenever it was possible, e.g. by 
comparing independent spectra taken at short time intervals be- 
fore merging them, or by comparing W values for those lines 
seen twice, on two consecutive orders. 

At McDonald Observatory echelle spectra of HD 34078 
were obtained with the Harlan J. Smith 2.7 m telescope. The 
strongest molecular features from CN, CH + , and CH are de- 
tected, as are the K line of Ca II and Ca I A4226, although the 
latter falls next to a CCD glitch and thus cannot be reliably ex- 
tracted. The data were reduced in the same fashion as before 
(cf. R03). A global multi-order fit of the entire CCD chip was 
performed for each Th-Ar spectrum, yielding residuals below 
0.001 A (or 0.08 km s _1 ). Measured radial velocities of (non- 
variable) absorption lines towards the comparison star £ Per 
show a scatter that is consistent with the residuals from the 
wavelength calibration. From the Th-Ar data we measured the 
instrumental resolution of our spectra, which turned out to be 
R = 170000. Stellar exposures were 30 m long for HD 34078, 
and typically 5 m for £ Per. 
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Fig. 5. Equivalent widths of CH A4300 measured in HD 34078 
(lower panel) and £ Per (upper panel) spectra (see Table |2). 
Open triangles correspond to raw OHP/ELODIE values, which 
appear systematically too low (filled triangles correspond to 
corrected values, see text). Tick marks numbered 1 to 8 show 
the epochs at which the FUSE spectra have been taken. 



A few observations of HD 34078 were also performed 
at other telescopes. Two spectra were obtained using the 
MAESTRO spectrograph, fed by the 2-m telescope at the 
Terskol Observatory (TE) in Northern Caucasus (the reso- 
lution was 120000). Three more spectra were obtained us- 
ing the fiber-fed echelle spectrograph installed at the 1.8-m 
telescope of the Bohyunsan Optical Astronomy Observatory 
(BOAO) in South Korea (some description can be found in 



Galazutd inov et al. 20051) . Modes providing resolutions, R = 



30 000 or 45 000 were employed (Table©). Finally, a spectrum 
was taken in service mode at the Cal ar Alto observatory with 
the FOCES spectrograph (R = 40 000; |Pfeiffer et alll998l) . For 
these additional runs, £ Per was not observed; we are neverthe- 
less confident that these data are homogeneous with respect to 
the other spectra obtained. 

3.2. Equivalent widths and column densities 
3.2.1 . Equivalent width estimates 

Equivalent width measurements were performed in a similar 
way as in R03 for CH A4300. Regarding W4 23 2(CH + ), the ac- 
curacy is limited by the poor definition of the continuum due to 
blending with a stellar line (cf. Fig. 7 in R03). Then, to avoid 
possible resolution-dependent effects on the estimate of W val- 
ues and to improve the homogeneity and accuracy of our mea- 
surements, we fitted this stellar line with a Gaussian. Its posi- 
tion (<5A = —0.364 A with respect to the CH + absorption) and 
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Fig. 6. Same as in Fig. but for the equivalent widths of 
CH+ A4232 in HD 34078 spectra. Raw OHP/ELODIE values 
(empty triangles) have been corrected by + 8% (filled triangles) 

FWHM (0.46 A) were fixed while the depth was varied so as to 
match the continuum blueward of the CH + absorption as well 
as possible (a typical value for the depth is 3% of the contin- 
uum). Once the spectra have been normalised in this way, we 
measure W and its uncertainty as done for CH A4300. This 
procedure was applied to all spectra, including those obtained 
prior to 2003; for the latter, this results in W values which dif- 
fer slightly from those given in R03. 

Some observations of HD 34078 were spread over 2 or 3 
consecutive nights (in particular the recent, high S/N, McD and 
OHP/SOPHIE spectra), which allowed us to check that no sig- 
nificant day-to-day variations are present; we then measured W 
on the co-added spectrum. For the same reason, one BOAO and 
one TE measurement obtained two days apart in March 2004 
were combined (thus Table|2]contains 25 entries). 

Uncertainties were estimated in a conservative way, in- 
cluding the error due to finite pixel-to-pixel S/N ratio and 
the one in continuum placement; we assume as in R03 that 
the two sources of errors combine quadratic ally. W esti- 
mates for HD 34078 and their associated errors are given in 
Table|2]for HD 34078 (CH A4300 and CH+ A4232) and ( Per 
(CH A4300); note that concerning OHP/ELODIE spectra, val- 
ues corrected as explained below are given for both HD 34078 
and ( Per. 

3.2.2. Consistency of all measurements 

In Fig. [5] (upper panel), we show results from the CH A4300 
( Per observations performed between 2003 and 2007, to 
whic h we add o lder measurements from ICrane et a l. (1995) 
and lAllenl (1994). These are consistent with a constant value 



of W43oo(CH). However, careful examination reveals a small 
offset of about -6% for the OHP/ELODIE values with respect 
to the other measurements. Since the same offset appears to 
be present in the HD 34078 W43oo(CH) values (lower panel, 
empty triangles), this effect is very likely due to scattered light 
in the ELODIE spectrograph (S. Ilovaisky and P. Prugnel, per- 
sonal communication). Thus, we applied a +6% correction 
(scattered light from the target does lead to a multiplicative cor- 
rection on W values) to all OHP ELODIE CH A4300 values 
(filled triangles), including those presented in R03. 

Unfortunately, the CH + A4232 line towards ( Per is too 
faint (W ~ 2.5 mA) to assess whether OHP/ELODIE mea- 
surements of this transition are also affected by scattered light. 
Then, to determine the correction for CH + A4232 (which may 
be different from that for CH A4300), we have to rely on the 
HD 34078 data themselves (Fig. [6}. By comparing the sets of 
OHP and McD values, we find that an offset of about +8% 
needs to be applied to the OHP/ELODIE values to bring both 
sets of points in good mutual agreement (filled triangles in 
Fig. 0. After correcting the OHP/ELODIE W values in this 
way, it is apparent in Fig. [5] and Fig. [6] that nearly simultane- 
ous measurements performed at different telescopes yield con- 
sistent values, within errors. This is a direct indication that un- 
certainties are not underestimated and that, after correction of 
the OHP/ELODIE W values, the whole set of data is homoge- 
neous. 

The OHP/ELODIE W measurements of CH+ A3957 are 
not accurate enough in comparison to those for CH A4300 
or CH + A4232 to be really useful (further, they are af- 
fected by scattered light in an unknown way). In contrast, 
the recent OHP/SOPHIE spectra provide good S/N values for 
W / 3957(CH + ). In Fig. [7] (panels c and d) we show the variation 
of W versus time for both CH + transitions since September 
2006; as can be seen, the two sets of measurements are very 
consistent. We also display for the same epochs the behav- 
ior of W 43 oo(CH) for both HD 34078 (panel a) and £ Per 
(b); the latter values remained constant while the variations 
seen for CH in HD 34078 are qualitatively similar to those 
observed for CH + . The whole set of values for W4232(CH + ) 
and W43oo(CH) shows a fairly smooth variation, with appar- 
ently little or no variations with timescales smaller than a few 
months. 

Fig. [5] Fig. [6] and Fig. [7] strongly suggest that the equiva- 
lent widths varied for HD 34078 while W43oo(CH) remained 
constant for £ Per. Let us now assess in a quantitative way the 
statistical significance of time changes in the observed W val- 
ues. To this end, we perform a \ 2 test on W / 43oo(CH) values 
for both the £ Per and HD 34078 sets of measurements, in or- 
der to check whether the assumption of a constant W can be 
accepted or rejected. Using the whole data set, we find that the 
C Per data are consistent with a constant value, W / 43oo(CH)= 
15.93 mA, while the assumption of a constant W value can be 
rejected at the 5.2 a level for HD 34078 (the W value which 
minimizes \ 2 is 52.3 mA). Note that these conclusions are un- 
changed if the OHP/ELODIE data are removed. Equivalently, 
if the distribution of x = (W— < W >)/cr values is com- 
pared to a Gaussian with a x — 1, one finds that it is consistent 
for £ Per values and inconsistent for HD 34078 (the x 2 test de- 
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Fig. 7. Equivalent width of CH A4300 (a), CH+ A4232 (c) and 
CH+ A3957 (d) in HD 34078 spectra since September 2006 
(symbols used have the same meaning as in Fig. [5j. The equiv- 
alent width of CH A4300 in £ Per spectra (b) is shown for com- 
parison (the horizontal line corresponds to the weighted aver- 
age of the displayed values). Note the similarity between both 
CH+ transitions and the good CH - CH+ correlation. 



scribed above is just a way to quantify these statements). We 
can therefore conclude that real time variations of W43oo(CH) 
occurred during the 2003-2008 period for HD 34078 (for CH+ 
towards HD 34078, the assumption of a constant W is rejected 
at the 3.2 a level). 

3.2.3. From equivalent widths to column densities 

Since the absorption lines considered here are not completely 
optically thin, we cannot infer N values directly from W ones. 
We will then rely on the highest resolution data, from McD 
observations. In the latter, the CH + A4232 profiles are fully 
resolved; one can thus derive the true optical depth and get 
iV(CH + ) by direct integration. Regarding CH, the determina- 
tion of TV and of the velocity profile is complicated by the 
struc t ure o f the ground level, due to A doubling. As shown by 
Lienl d 1984ft . neglecting this effect may cause an underestimate 
in A^(CH) and a broadening of the profile, when the intrinsic 
width is small enough (i.e. comparable or smaller than the split- 
ting of the ground level which corresponds to 1.43 km s _1 ). 

CH profiles include at least two components (a strong 
asymmetric and narrow one superimposed onto weak shallow 
absorption) and fitting with Voigt profiles remains somewhat 
arbitrary (Fig. 5 and Fig. 7 in R03). Thus, to get the intrinsic 
CH pixel optical depth profile without any a priori decomposi- 
tion and perform the detailed comparison with the CH + profile 
allowed by the quality of the McD data, we use a Bayesian 



inversion procedure as done bv lPichon et al. I (l200ll) in the con- 
text of Lya absorption in QSO spectra. Equally populated sub- 
levels are assumed (we checked that the relative strength of the 
CH A3886 and CH A3890 lines in the O HP/SOPHIE spectra is 
consistent with this hypothesis; see also Lienl 1984 , for other 
lines of sight). We find in the end that taking into account A 
doubling induces corrections on A^(CH) which are no larger 
than 1.3%. 

Significant changes in the CH and CH + profiles are seen 
(see Sect. 13.3b but these remain relatively small and further, 
the optical thickness of either CH A4300 or CH+ A4232 does 
not exceed 1 . One can therefore expect a simple empirical rela- 
tion to hold between W and N, allowing us to infer N values 
from W measurements performed from lower resolution data. 
In Fig. [8] we show W/X versus NfX for both CH A4300 and 
CH + A4232 from McD spectra. As can be seen, both sets of 
points are well fitted by a single straight line (dotted line), as a 
result of the similarity of CH and CH + velocity profiles (cf R03 
and Sect. 13.31 . This fit corresponds to the following relations, 



(1) 

(2) 



JV(CH) = 2.42 x 10 n W 1AS 
and 

AT(CH+) = 2.29 x lO 11 ^ 1 - 48 



for CH A4300 and CH+ A4232 respectively, with W in mA 
and N in cm -2 . It is noteworthy that although we made no as- 
sumption on the shape of line profiles in the above procedure, 
the best fit turns out to be close to the curve of growth for a 
single gaussian component with b = 3.2 km s . All measure- 
ments are nicely bracketted by curves of growth with b values 
between 2.8 and 3.5 km s . We next assume that the N(W) 
best fit relation holds at all epochs and use it to infer N from W 
for measurements other than the McD ones. In Table|2] we list 
A^(CH) and N{CW + ) values; for McD spectra, these are drawn 
directly from line profiles while for lower resolution data, N 
is obtained from W through Eqs.[T]and[2] The resulting uncer- 
tainty in N should in principle involve two terms, one related 
to the errors on W and the other to the scatter about the best fit; 
in practice, the latter appears to be small and we ignore it. 

3.2.4. CH and CH+ column density variations 

Fig. |9] shows the variation of Af(CH) and 7V(CH+) observed 
between 1998 and 2007. In R03 we found that W4 30 o(CH) 
increased from 1990 to 2002; the +6% applied to the 
OHP/ELODIE data results in a larger amplitude for this varia- 
tion (21% instead of 14%). It now appears that A^(CH) reached 
a maximum during the 2000-2002 interval, and has been de- 
creasing since then to reach in 2006 a value similar to those 
observed before 1998. To better characterise the long-term 
variation, we smoothed the data using a gaussian window 
with a FWHM of one year. In the averaging, each N value 
is weighted according to l/a 2 (N) where cr(N) is the un- 
certainty; the resulting curve (together with ±3<t bounds) is 
shown in Fig. [9] The typical timescale for these variations, 
tjv, defined as the time needed for A" to change by 10% (i.e. 
t n = [dN/iOANdt)}- 1 ) is about 1.5-2 yr. 
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Fig. 8. Equivalent width versus column density for high reso- 
lution observations (McDonald). CH (filled squares) and CH + 
(empty squares) absorption lines span different ranges in opac- 
ity which helps to constrain the best linear fit (dotted line: 
y = -15.232 + 0.676 x). A curve of growth with b = 3.2 
km s _1 also provides a good fit to the data (plain line; the two 
additional lines correspond to b = 2.8 and b = 3.5 km s _1 and 
bracket the data points well). The large b curve (i.e. optically 
thin limit) is shown (dashed line) 

The new data clearly indicate that additional variations, 
with shorter timescales, are also present as was suspected from 
earlier data (R03). This is especially clear in the most re- 
cent 2006-2008 results: N(CH) decreased by -10% between 
September 2006 and February 2007, (tjv — 0.5 yr) and then in- 
creased by 9% to reach the last February 2008 value (see Fig. [7] 
for a zoom on the recent W data). 

We now consider the CH + measurements shown as 
empty squares in Fig. [9] The +8% correction applied to the 
OHP/ELODIE data somewhat affects the previous conclusion 
of R03 of a possible decrease of 7V(CH + ) between 1990 and 
2002: 7V(CH + ) now looks essentially constant over the long 
term, with no large amplitude variations as those seen for CH. 
Beyond 2000, the general pattern is relatively similar for both 
species. However, although the long-term variations in CH and 
CH + seem to be loosely correlated, there is good correspon- 
dence between the short-term ones; this is especially clear in 
Fig. [7] (note in particular the coincidence of the local min- 
ima in February 2007 for both species). Although the varia- 
tion patterns are qualitatively similar, the amplitudes are sig- 
nificantly different. For instance, the increase of iV(CH + ) be- 
tween February 2007 and 2008 (+26%) is notably larger than 
that of iV(CH) (+9%). 

We shall now focus on the high resolution spectra of the 
CH A4300 and CH+ A4232 absorption lines in order to investi- 




2000 2005 
Time (year) 



Fig. 9. Evolution of the CH and CH + column densities with 
time. The solid thick lines show the long-term evolution ob- 
tained after smoothing using a gaussian window with a FWHM 
of one year (dashed lines indicate ±3er bounds on smoothed 
values). 

gate whether the observed long and short-term variations in W 
are accompanied by detectable changes in line profiles. Profile 
variations might also be present which are not necessarily re- 
flected in appreciable W changes. 

3.3. High resolution CH and CH + line profiles and 
their variations 

The high resolution profiles obtained at McDonald Observatory 
have all been corrected for the Doppler shift due to the Earth 
velocity and have been brought in the LSR system. Yet, after 
correction, the wavelength of the £ Per CN and CH lines still 
show slight fluctuations in position (of at most a few mA) from 
one epoch to another (given the excellent S/N and the sharp 
line profiles, misalignments by only a few mA are sufficient 
to induce significant profile differences). These fluctuations are 
identical for both CN and CH lines to within 1 mA and show a 
good correlation with those of the HD 34078 lines. They must 
then be due to some inaccuracy in the LSR correction from 
one epoch to another. To improve the accuracy in the align- 
ment of the HD 34078 absorption lines, we assume that ( Per 
lines have been stable in position (this is confirmed by SOPHIE 
spectra whose wavelength scale is accurately defined) and infer 
the value of the relative shifts for each epoch. These are used 
to slightly adjust the position of HD 34078 lines. 

The CH and CH + line profiles are quite similar (Fig.fTUb: 
both include a narrow component with a full width at half max- 
imum of about 5.5 km s _1 (corresponding to a b parameter of 
3.3 km s" 1 ) and shallow absorption extending from —4 up to 
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Fig. 10. Comparison of McD CH A4300 (lower left) and 
CH + A4232 (upper left) profiles observed in January 2003 and 
October 2005. During this 2.81 yr time interval W has varied 
by -11.7% and +0.3% respectively. Same for OHP/SOPHIE 
CH A4300 (lower right) and CH+ A4232 (upper right) profiles 
obtained in September 2006 and February 2007 (W variation 
of -7.1% and - 6.6 % over 0.43 yr). Note the good stability of 
the red side of the profiles. 



17 km s _1 . The CH and CH + narrow components cover the 
same velocity range (~ 2 - 11 km s _1 ) but their shapes are 
significantly different: CH displays a steeper blue edge, while 
the opposite holds for CH + . This results in a shift of about 
+1.6 km s" 1 for the CH + line centre with respect to that of 
CH. It is noteworthy that the velocity range covered by the nar- 
row CH or CH + absorption coincides quite well with the range 
over which CO emission is observed. 

Since we are likely probing molecular gas closely associ- 
ated with HD 34078 (i.e. with peculiar physical conditions), 
one may wonder whether the CH and CH + profiles show sig- 
nificant differences with those seen on other lines of sight. We 
simp ly note that the presence of a we ak broad component is 



rare (ICrane et al.lll995t ICrawfordl ll995) but such shallow ab 



sorption would be difficult to detect in most spectra; one of 
very low amplitude (~ 1% instead of ~ 5% in our spectra) is 
however present in the very high S/N spectrum of ( Oph pre- 
sented by lCrane et"aD dl99lh . 
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Fig. 11. Comparison of the HD 34078 OHP spectra ob- 
tained in February 2007 (black) and February 2008 (red) for 
the CH A4300 (upper left), CH+ A4232 (lower right) and 
CH + A3957 (upper right) lines. Note the good consistency be- 
tween the two CH + transitions. During this 1.0 yr interval, the 
( Per CH A4300 profile remained stable (lower left). 



By comparing successive spectra, we find that fluctuations 
in W are most often associated with changes in the blue and 
central part of the profiles, the red side suffering very little 
variation. This is apparent in Fig. [To] (left panels) in which we 
display CH McD profiles for the January 2003 and October 
2005 epochs. During this 2.8 yr time interval W43oo(CH) has 
decreased by 12% and the FWHM of the narrow component 
decreased by 5%. W4232(CH + ) remained nearly constant over 
the same period but significant profiles changes are neverthe- 
less clearly present; the line is slightly deeper, which com- 
pensates for a decrease in FWHM comparable to that seen 
for CH. Another example is shown in Fig. [10] (right panels), 
where we compare OHP/SOPHIE CH and CH + profiles taken 
over a much shorter interval (0.43 yr from September 2006 to 
February 2007). In this case, W 43 oo(CH) and W 4 232(CH + ) have 
decreased by a comparable amount (~ 7%) while both profiles 
became slightly narrower. Thus correlated W variations ap- 
pear to be associated with similar profile changes. Since 2003, 
W values show no systematic trend but rather erratic fluctua- 
tions; the same is true for profiles and generally, for two epochs 
with comparable W values, the profiles are quite similar (the 
CH+ A4232 January 2003 and October 2005 profiles being an 
exception). In Fig.QT]we display the recent February 2007 and 
2008 CH and CH + profiles; the latter show a marked increase 
in W (cf Fig.|7), corresponding to +26% for N(CH+) t (Note 
that during the same interval, the CH A4300 ( Per profile re- 
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Fig. 12. Comparison of the FUSE spectra obtained in January 
2000 (LWRS, black line) and October 2004 (MDRS, red line) 
in the region around 1050 A (v and J values for the lower level 
of each H2 transition are indicated; the HD line is from J = 0). 
In the MDRS spectrum, narrow lines tend to appear deeper and 
the damped (0,0) and (0,1) features slightly broader. 

mained stable.) We do not see appreciable variations of the 
broad shallow component, but given its weakness, the signif- 
icance of this result is limited. 

Unfortunately, the S/N ratio for the CN line profiles of 
HD 34078 is not sufficient to allow a search for variations with 
a sensitivity comparable to the one attained for CH and CH + . 
Further the strongest CN line is clearly affected by variations 
of blended stellar (C IV) absorption. Thus we shall not discuss 
variations of CN features. 

4. Variations in the H 2 column density 

Bois se et"al] (2005) analyzed the first 5 FUSE spectra taken 



since 2000 and detected no variation for the H2 column den- 
sity with an upper limit of 5% at a 3a confidence level. Three 
additional spectra were obtained by FUSE in September 2003, 
February 2004 and October 2004 which allow us to follow the 
time behavior of H2 absorption during nearly five years since 
January 2000 and then probe the spatial distribution of the gas 
over scales up to 104 AU. 

4. 1. Contamination of HD 34078 spectra by nebular 
emission 

Comparing the October 2004 (8 th ) spectrum to the previous 
ones, we find that it shows noticeable differences (Fig.fTZl): nar- 
row lines are deeper and damped H 2 lines are slightly broader 
[as if N(H 2 ) had increased]. This brought to our attention a 
potential difficulty that had not been considered in B05: con- 
tamination of the HD 34078 spectrum by diffuse emission 
from the IC 405 nebula. Indeed, the 8 th spectrum was obtained 
with the MDRS aperture while all others were taken using the 
larger LWRS aperture. Given the difference in size (MDRS: 4 



x 20 arcsec 2 ; LWRS: 30 x 30 arcsec 2 ) and the intense dif- 
fuse emission detec ted close to the HD 34078 line of sight by 
France et alJ ((2004), the peculiarities of the 8 th spectrum might 
just be due to a lower level of contamination of the spectrum 
by diffuse emission. In Appendix A, we estimate the nebular 
contribution to LWRS spectra and conclude that it can explain 
the difference between FUSE spectra 1 to 7 and the 8 th MDRS 
one. Thus only spectra 1 to 7 will be considered below in our 
search for variations in H 2 lines. 

The importance of diffuse emission contamination in 
LWRS spectra also implies some limitation in our search for 
variations: the aperture may not be located exactly at the same 
position on the sky at all epochs resulting in a slightly vari- 
able contribution from diffuse emission if gradients are present. 
Note that since the diffuse to stellar flux ratio decreases with 
wavelength, H 2 systems at longer wavelengths are best suited 
to minimize the contamination by diffuse emission. Regarding 
the study of the gas properties towards HD34078, the 8 th spec- 
trum is clearly to be preferred for two reasons: i) it should be 
much less affected by diffuse emission and ii) the S/N ratio is 
significantly higher than for previous spectra due to an integra- 
tion time (22 500 s) about four times longer than at epochs 1 
to 7 (~ 6000 s). A redetermination of the gas properties based 
on the 8 th spectrum (H 2 excitation diagram in particular) will 
be presented elsewhere; we simply note here that the detection 
of absorption lines from all excited H 2 levels quoted in B05 is 
confirmed. 



4.2. Variations in N(H 2 , J=0) 

As in B05, we perform a direct comparison of the LWRS spec- 
tra, after relative flux intercalibration and adjusment of the 
wavelength scale. This procedure is applied independently to 
three portions of the spectra located at about 1050, 1063 and 
1078 A, corresponding to the (4-0), (3-0), (2-0) H 2 Lyman 
bands respectively. Each of these broad features is a blend of 
four H 2 lines arising from the J = 0, 1 and 2 levels. A good 
relative flux calibration is easily obtained (as for spectra 1 to 5), 
indicating that the shape of the stellar spectrum does not vary 
(known stellar li nes for such Q9 star s are indeed weak and rare 
in these regions: iPellerin et al ] l2002h . Using narrow high J H 2 
lines adjacent to the broad H 2 absorptions of interest, we get an 
accuracy in the wavelength alignment of about 0.01 A for the 
1050, 1063 and 1078 A absorption systems. 

We now focus on the blue edge of each broad H 2 system 
which presents good sensitivity to changes in Af(H 2 , J = 0). 
A zoom of this region for the 1050 A system is shown in 
Fig. [T3l (upper panel). All spectra, corrected in flux and wave- 
length as described above, are superimposed. They are all sim- 
ilar and an average spectrum can therefore be computed (thick 
line). The difference (Ali) between one individual spectrum, 
i, and the mean is displayed in the lower panel for each epoch 
(i = 1, 7 from top to bottom). The 3<r dispersion on Al, 
among the seven epochs is indicated as a function of wave- 
length (dashed lines). Away from the J = line, the AI profile 
is consistent with no variation. In the region close to the J = 
line (displayed in red) where variations in N( J = 0) would 
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Fig. 13. FUSE LWRS HD 34078 spectra near the 1050 A broad 
(4-0) H2 Lyman band. Upper panel : all epochs, after flux in- 
tercalibration and alignment (see text). The mean spectrum is 
shown (thick line); tick marks indicate the features shown in 
Fig-ERtwo are from H2, the third is from Ar I). Lower panel : 
difference between the spectrum for each epoch and the mean 
spectrum (epochs 1 to 7, shifted from top to bottom for clarity; 
same scale as in the upper panel). The dashed line shows the 
3cr variance among the different epochs. The region used in the 
analysis is shown in red. 

induce changes in the profiles, spectra 1 to 7 are also con- 
sistent with the mean spectrum. Similarly, B05 have adjusted 
the first spectra using log N( J — 0) = 20.52 and concluded 
that the variation among the five first spectra was lower than 
5%. Indeed, an increase (decrease) of this amplitude roughly 
corresponds to a difference that follows the lower (upper) 3<r 
profile in Fig. 10. We conclude from our analysis, that /V(H2) 
changed by less than 5% at the 3a level between January 2000 
and February 2004 while A^(CH) has undergone variations as 
large as 20% over the same time interval (cf Fig. 6). 

5. Observations summary 

5. 1. The stability of( Per CH absorption 

As mentioned above, £ Per was observed primarily to test 
the homogeneity of our measurements. After correction of the 
small offset found in the OHP/ELODIE data, the whole set of 
W values appears remarkably consistent. The latter, and the 
most recent (2006 - 2008) data in particular, lead us to a strong 
twofold conclusion: i) W43oo(CH) remains constant towards 
C Per and ii) N(CH) and N(CH+) do vary towards HD 34078. 
Indeed, if variations of instrumental origin were responsible for 
changes in HD 34078 values, one should invoke a very unlikely 
"conspiracy" to explain the stability of £ Per lines. The fact that 
the same behavior is observed for distinct CH + lines is a strong 
additional proof of the reality of HD 34078's variations. 

£ Per was observed during 5 years and over this time inter- 
val, the drift of the line of sight through the foreground cloud 
has be en significant. The dist ance to the cloud is thought to be 
350pc dHilton & Lahullalll995|) and thus the drift of the line of 



sight amounts to 17.8 AU. The constancy of N(CH) then shows 
that over this scale and below, there is no marked structure in 
the cloud. The 3a upper limit on relative variations of AT(CH) 
is about 6%. We derive this value simply from the raw aver- 
age and rms scatter of McD and OHP/SOPHIE measurements, 
assuming the CH A4300 line is optically thin; a more detailed 
analysis of the ( Per data, including CN and CH + lines, will 
be presented elsewhere. 

5.2. Main observational results on HD 34078 

From multiwavelength observations of the gas towards 
HD 34078 and in its close environment, we get the following 
results: 

★ the 12 CO(2-l) emission map of the HD 34078 field shows a 
pronounced peak coincident with the star's position, clearly 
indicating that molecular gas seen in absorption is closely 
associated with HD 34078. The extent and morphology of 
the CO em ission correla t es wel l with the 24 pm dust emis- 
sion arc of Fran ce et al. (2007). Presumably as a result of 
the interaction between HD 34078 and the ambient cloud, 
the CO(2-1)/CO(1-0) ratio is anomalously large, pointing 
towards dense (10 3 — 10 4 cm' 3 ) and warm (T > 12K) 
emitting gas and further, a remarkable kinematical pattern 
with doubled-peaked profiles is observed, 

★ we confirm the reality of rapid, large amplitude (typi- 
cally 10% yr _1 ) and correlated variations of A^(CH) and 
iV(CH + ) towards HD 34078. The velocity ranges covered 
by CH and CH + narrow absorption coincide well with that 
of CO emission. Variations in CH and CH + line profiles 
are unambiguously detected; these occur mainly in the blue 
part of the narrow absorption. A broad shallow and rela- 
tively stable component is seen for both CH and CH + in 
the interval [-4, 17 km s _1 ], 

★ comparison of LWRS and MDRS FUSE spectra reveals 
that the 7 LWRS spectra available are significantly contam- 
inated by diffuse light from IC 405. The absence of varia- 
tions in the LWRS profiles of H2 J = lines yields a 3a 
upper limit of 5% on N values, extending the result of B05 
over nearly four years (or 90 AU). 

6. Discussion 

Given the marked contrast between the stability of ( Per CH 
lines and the rapid, large amplitude variations seen for CH and 
CH + towards HD 34078, we shall assume in the following dis- 
cussion that these variations can be attributed entirely to phe- 
nomena associated with the star/cloud interaction and not to 
small scale structure in cold gas. 

6. 1. Towards a coherent picture of the close 
environment of HD34078 

From the broad set of observations available, a coherent sce- 
nario emerges which can be summarized as follows. HD 34078 
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rec ently encounte red a molecular cloud, as originally suggested 
by iHerbigl dl958f) . The stellar wind impacts the ambient ma- 
terial resulting in a shell of compressed, highly excited gas. 
Modelling of both the H 2 excitation (B05) and CO emission 
(Sect. 2.2) provides consistent estimates for the density in the 
shell, n ~ 10 4 cm -3 . Dust grains located in this region are 
directly exposed to the intense UV flux from the O star and 
strongly emit at infr ared wavelengths. T hus, both the mid-IR 
emission detected by lFrance et al.l (120071) and the CO emission 
probably delineate the part of the shell seen edge-on by the ob- 
server, accounting for the good correlation between the Spitzer 
8 or 24 /xm arc and the CO map (Fig. 1). 

B05 favored a two-cloud model on the basis of the di- 
chotomy in the physical conditions derived for the highly ex- 
cited H 2 on the one hand and for all other molecular absorp- 
tion on the other hand. They discarded a single cloud model 
(although it would have naturally explained the similar ve- 
locity for all absorption lines and the presence of preexist- 
ing diffuse H 2 near HD 34078 ...) because they implicitly as- 
sumed that molecules should be photodissociated at the very 
small distance (a few 0.01 pc) implied by the modelling of 
the H 2 excitation. In fact, this argument is valid only in a 
model that is stationary regarding the formation/destruction of 
molecules. Given the large space velocity (V s ) of HD 34078, its 
arrival is so recent that probably, no such steady-state equilib- 
rium could be established. Rather, as the O star is approaching 
the cloud, a photoionisation and photodissociation front devel- 
ops, moving at velocitie s of the order of a few km s _1 only 
( Bertoldi & Drainall996 ). i.e. well below the star velocity. The 
distance between HD 34078 and these fronts then gradually 
decreases and the velocity of the latter gets higher (the front 
velocities increase with stellar flux), up to a point where the 
star and front velocities equate. 

Moreover, the stellar wind has a strong mechanical im- 
pact on the surrounding gas and the latter is gradually set into 
motion as a result of momentum flux. When the star is close 
enough to the cloud, a stationary bow shock is established at 
a position nearly coincident with that of the photodissociation 
front. At this stage, the distance (i?o = AS, see Fig. \14\ be- 
tween the star and the apex of the shock is well defined and 
remains constant as far as the density of the ambient cloud and 
the wind properties (mass loss rate and terminal velocity) do 
not change. In this picture (Fig. [Pil l, whether the dynamical 
steady-state regime is established or not, the molecular material 
located beyond the front/shock surface should be little affected 
by the presence of the closeby star, thereby accounting for the 
characteristics of the molecular components that B05 assigned 
to the "translucent" component. 

We conclude that a model involving a single cloud located 
very close to HD 34078 may, in fact, be consistent with all ob- 
servations. A key issue remains open however: at which stage 
of the star-cloud interaction have we captured HD 34078 and 
its close environment ? A stationary dynamical regime is ex- 
pected to establish ultimately and in the thin-shell limit, de- 
tailed models describing the geometry and v elocity structure 



Ambient cloud 




Fig. 14. Cartoon showing some geometrical parameters associ- 
ated with the bow shock around HD 34078 (h is the apparent 
standoff distance; h ~ 0.04 pc). The velocity vector V s shown 
as the bold arrow is inclined by 8 — 30° with respect to the 
plane of the sky. "S" marks the star position while the dot- 
dashed horizontal line represents the line of sight to the star 
(the observer being on the left). The double curve is the pro- 
jection of the proposed bow shock. Optical and UV spectra of 
HD 34078 probe the bow shock/hot PDR at L 2 (with an incli- 
nation angle ip\ we estimate SL2 — 0.1 pc) and the foreground 
ambient cloud while IRAM 230 GHz observations involve all 
material enclosed into the beam. 



and observations (Appendix B). From this analysis, we con- 
clude that the IR/CO arc does not display the properties ex- 
pected for a stationary bow shock and that we are possibly ob- 
serving a "nascent bow shock", i.e. the wind/cloud interaction 
at an early evolutionary stage, well before the formation of a 
steady-state flow around the star. We also find that the geomet- 
rical constraints provided by the IR data, h ~ Rq ~ 0.04 pc, 
are roughly consistent with the radiation field implied by the 
modelling of H2 excitation (B05). 



A variant of the above scenario is suggested by lFrance et al 



of steady-state b ow shocks are avai lable (e.g. Ivan Buren et al 



1990tlMac Low et al][l99llWilkinlll996h . We thus performed 
a detailed comparison between the predictions of these models 



(2004), who proposed that differential extinction is present be- 
tween HD 34078 and the surrounding diffuse emission, in or- 
der to explain the increase of the diff use to stellar ratio a t far- 
UV wavelengths. More specifically, France et al.l (120041) pro- 
pose that a small clump lies in front of HD 34078 (Fig. [T5T >. 
With an extent no larger that about 20 ", the latter could in- 
duce the observed HD 34078 extinction without affecting the 
surrounding nebular emission. 

Such a picture is attractive in our context because 20 " is 
approximately the size of the area over which a "dip" is seen in 
the CO profiles, suggesting that the double-peaked line shapes 
might be due to narrow absorption rather than to velocity struc- 
ture in the emitting gas. In this scenario, CO emission could 
originate from the ambient cloud background to HD 34078 - 
especially its outer boundary, compressed by the wind - ac- 
counting for the widespread emission seen over most of the 
field (Fig.rT3Tl. Additional emission could come from the clump 
itself, in particular from the region located immediately beyond 
the hot PDR facing the star, explaining the enhanced emission 
close to HD 34078 and the high excitation of the emitting ma- 
terial. Gas located on the cold side of the clump facing the 
observer, should be little affected by the interaction, as in the 
bow-shock scenario. This region corresponds to the translucent 
component in B05's model; with a small CO excitation temper- 
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Fig. 15. Alternative scenario to the one presented in Fig. [14] 
Here, HD 34078's wind slightly modifies the shape of the am- 
bient cloud boundary, giving rise to the IR arc. The clump red- 
dens the star but not the surrounding nebular emission; further, 
gas located on its cold side (facing the observer) induces nar- 
row absorption in the 12 CO(2 — 1) emission profiles as well as 
some visible and UV absorption lines. Excited H 2 lines and a 
significant fraction of CH and CH + absorption originates from 
the PDR located on the hot side of the clump, intersected by 
the line of sight at L 2 . 



ature and low velocity dispersion, this gas could induce narrow 
absorption in the background CO emission. The close similar- 
ity of the CO dip velocity at the star position, V ~ 6.0 km s _1 
(this value is stable to within 0.2 km s _1 among spectra dis- 
playing the dip), and the velocity of maximum narrow CH and 
CH + absorption (see Fig. ITOb is consistent with absorption be- 
ing responsible for the dip in CO profiles. The clump should be 
located at a distance from the star of about 0. 1 pc (comparable 
to that of point L2 in the bow-shock picture) so that the hot 
PDR facing the star is exposed to a radiation field large enough 
to account for the H2 excitation (Appendix B, B05). In fact, 
the bow-shock and clump scenario are more or less equivalent, 
both of them involving the presence along the line of sight of 
a shell of dense gas illuminated by intense UV radiation; the 
main difference is that in the clump picture, the ambient mate- 
rial is not distributed in a continuous manner around HD 34078. 

6.2. CH and CH+ abundance and time variations 

In the picture described above (Fig. [14), the dense shell at 
the wind/cloud interface is nearly static in the ambient cloud's 
frame and as time elapses, the point (L2) at which the line 
of sight intersects the shell drifts over the latter towards the 
apex {A) at a velocity Vj[/ cos tp, where V[ is the transverse 
velocity and <p is the inclination of the shell at L 2 - Similarly, 
in the clump scenario (Fig. |T3T >. the line of sight drifts over 
the wind/clump interface. The rapid CH and CH + variations 
observed imply that a significant fraction of these species is 
enclosed in a very localised region (with a size of about 10 
AU, corresponding to a time interval of 6 months) and in our 
context, it is natural to assume that time variations arise from 
structure over the shell of compressed gas intersected by the 



line of sight. Further, as already outlined by B05, the CH/H2 
ratio towards HD 34078 is anomalously large: for A^(CH) ~ 
10 14 cm -2 (the largest values reached, in 2000 - 2003) we get 
a ratio of 1.6 x 10 -7 , which appears to be 3 .7 times larger than 
the value inferred from the best fit given by lSheffer et alJ Q2008) 
(a similar result is obtained by com paring HD 34078 values to 
those compiled bv lWeltv et al . 2006). Since the CH - H 2 corre- 
lation is quite good, such a de viation is very signifi cant and in 



deed, in the Fig. 8 presented bv lSheffer et aD(l2008l) . HD 34078 



is clearly an outlier (note that surprisingly, HD 37903, the other 
star showing highly excited H2, has an anomalously low CH/H2 
ratio !). The CH + /H 2 ratio displays much more scatter among 
all sightlines, but nevertheless, N(CH + ) towards HD 34078 is 
amon g the largest values for Af(H 2 ) ~ 6 x 10 20 cm~ 2 (Fig. 



10 in lSheffer et all 120081) . In the end, the CH/CH+ ratio in 
HD 34078 is well within the rang e observed towards other stars 
(cf Fig. 10 in [Weltv et al. 2006). We thus apparently have a 



comparable relative excess of both CH and CH + , and again, 
it is natural to assume that the overproduction of these species 
occurs in the compressed shell. To induce such a large devi- 
ation in the CH/H2 ratio, the overproduction of CH must be 
quite large since molecular gas located beyond the shell prob- 
ably displays a more standard ratio. Similarly, the spatial dis- 
tribution of CH and CH + over the shell must be very inhomo- 
geneous at scales of about 10 AU to induce the observed varia- 
tions. Since the CH and CH + variations are strongly correlated, 
a single mechanism must be at work to explain the overproduc- 
tion of both species. 

Clearly, the CH and CH + time variations cannot be at- 
tributed to pure density structure. Indeed, no corresponding 
changes have been seen for iV(H 2 ) and further, explaining CH 
or CH + variations by more or less spherical clumps with a 
size of about 10 AU would require very high volumic densi- 
ties, n(H2) ~ 10 5-6 cm~ 3 , as argued by R03. In such gas, 
CH + would be rapidly destroyed by reactions with H 2 . Then, 
the structure is likely to be more chemical in nature. 

The production of CH + , which is not expected at thermal 
equilibrium, together with the correlated variations of CH and 
CH + , suggest that both molecules are mainly formed in the 
dense shell, through a MHD shock where the drift velocity be- 
tween ions and neutrals trigger the formation of both CH + and 



CH (jPineau des Forets et alJl986tlFlower & Pineau des Forets 
119981) . CH + could also be overproduced at the interface be- 
tween the ambient cloud and warmer gas, as suggested by 
I 1 J — l I 1 \ — L 

Duleyetal. (1992), Crawford (1995) and more recently by 
iLesaffre et al.1 (120071) . In the specific conditions prevailing 
around HD 34078, one may imagine that Rayleigh-Taylor in- 
stabilities develop efficiently (e.g. as a result of fluctuations in 
the wind properties or in the ambient cloud), leading to the for- 
mation of a turbulent mixing layer wit h pronounced small-scal e 
structure. HD 34078 being variable dMarchenko et al.1 1 19981) . 
stellar flux variations can also trigger or amplify the formation 
of small-scale structure at the interface and in the PDR. In re- 
gions hot enough to form significant amounts of CH + via the 
endoenergetic C + + H2 reaction (whatever the heating mech- 
anism, shocks or turbulent dissipation), small spatial temper- 
ature fluctuations will result in appreciable variations in the 
CH+ formation rate and then in the relative abundance of CH + 
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and CH (recall that CH is easily formed from CH + once the lat- 
ter species is present). One can indeed estimate that since the 
CH + formation rate scales as exp(— 4640/T), a local variation 
as small as 22 K around T ~ 1000 K is sufficient to induce 
a fluctuation of 10% in the local abundance. The broad, shal- 
low CH and CH + absorption components (Fig.fTOb could just 
be a signature of the highly turbulent velocity field at the inter- 
face (note the excellent agreement in the velocity intervals [-4 
to 17 km s _1 ] covered by the CH and CH + broad components, 
indicating that the species responsible for that component are 
cospatial). 

CH and CH + line profiles vary mainly on their blue side, 
at Vlsr — 4 — 5 km s~ 1 (Fig.[T0l>, suggesting this velocity for 
the gas lying at the interface near This material is there- 
fore blue-shifted with respect to the ambient cloud (at Vlsr 
~ 6.5 km s _1 ); this is consistent with the stellar wind push- 
ing ambient molecular gas towards the observer. The higher 
excitation derived for the CO component at Vlsr — 5 km s _1 
also nicely fits this view. Thus, the scenario described above 
might, at the same time, explain in a coherent way the large 
abundances of CH and CH + as well as their time variations. 

The amounts of excited H2 (J = 3 to 5) and CH + are 
known to correlate dLambert & D anks Il986h and is taken as 
an indication that the same mechanism is responsible for these 
species, the production of which requires an input of additional 
energy of yet unknown nature. Towards HD 34078, a marked 
excess of both excited H2 and CH + is observed, in agreement 
with the correlation seen among more standard lines of sight. 
We simply note that in our scenario, the large amount of ex- 
cited H2 and CH + is a direct consequence of the proximity of 
HD 34078, through its wind and UV flux. 

6.3. Small scale structure in quiescent H 2 gas 

In the picture that we propose, H2 gas located beyond the 
photodissociation front and shocked region is essentially unaf- 
fected by the presence of the star (dust grains should be some- 
what warmer due to the proximity of the star, but at densities of 
a few 10 2 cm -3 , this has little impact on the gas temperature). 
This part of the cloud gives by far the dominant contribution to 
the H2 column density in the J = level (in B05's model, the 
hot PDR represents less than 1 % of the total N value); it should 
also contain a significant fraction of the CO responsible for the 
UV absorption, about 1/4 of CH molecules (cf above) and most 
of the CN. Thus, the lack of variation in N(U 2 , J — 0) at a 
level better than 5% implies that no marked small-scale den- 
sity structure is present within the fraction of the ambient cloud 
probed by the drift of the line of sight between 2000 and 2004. 
If this region is representative of quiescent diffuse molecular 
material in general, the structure seen elsewhere for other trac- 
ers like H2CO, HCO + , and OH would be mainly "chemical" 
structure, possibly reflecting the specific formation/destruction 
processes relevant to these species. 

We further note that the lack of structure in quiescent H2 
gas is consistent with the stability of CH absorption lines to- 
wards £ Per, provided the CH/H2 abundance ratio is uniform 
at scales of about 10 AU within this quiescent cloud. 



7. Conclusions and prospects 

• By mapping 12 CO(2 — 1) emission around HD 34078, we 
have unambiguously shown that the molecular material seen 
in the foreground is c losely associa ted with the star, support- 
ing the suggestion by lHerbigl dl958h that HD 34078 currently 
encounters a molecular cloud. Repeated CH and CH + observa- 
tions, performed using the star £ Per as reference, confirm the 
reality of rapid and large amplitude variations of 7V(CH) and 
iV(CH+) along the line of sight. 

• The results altogether strongly suggest that the recent ar- 
rival of HD 34078 near the southern edge of a molecular cloud 
has given rise to a shell of dense gas at the interface between the 
stellar wind and H2 gas, the latter material belonging either to 
the distorted boundary of the a mbient cloud or to a small fore- 
ground clump, as suggested by France et alJ (120041) . The loca- 
tion of this shell relative to the star is consistent with constraints 
derived from earlier modelling of H2 excitation. By compar- 
ing the geometrica l characteristics of the IR arc detected by 
Franc e et al ] d2007h and the velocity field inferred from our 
CO or optical observations to predictions of steady-state bow 
shock models, we find that the latter are inconsistent with the 
observed properties. Therefore, we may be seeing this region at 
an early phase of the wind/cloud interaction, with a dense layer 
formed at the interface but no stationary flow yet established. 

• We propose that the large relative CH and CH + abun- 
dances originate from significant overproduction of CH + in the 
dense shell, due to the presence of a strong C-shock and/or of 
mixing of warm ionised gas and highly excited molecular ma- 
terial at the wind/cloud interface. The pronounced, correlated 
CH and CH + variations would then reflect marked chemical 
structure in the dense shell, possibly resulting from instabili- 
ties occuring at the interface. 

• No variations of N(H2, J = 0) have been found at a 
level of 5% (3er limit), extending the result obtained in B05 to 
a time interval of 4 years, or 110 AU. In the scenario that we 
propose, J = H2 molecules are mainly located beyond the 
dense shell. This indicates that beyond the photodissociation 
and photoionisation fronts, where the molecular material is 
not yet affected by the interaction with HD 34078, no marked 
small-scale scale structure is present and that the bulk of the 
mass is distributed relatively uniformly within the cloud. 

Let us now discuss some prospects about the time evo- 
lution of HD 34078's environment, assuming the scenario 
sketched in this paper is roughly correct. One may in partic- 
ular wonder whether the future evolution can induce signif- 
icant observable changes in the coming years. Using the es- 
timates for the observed standoff distance and its expected 
steady-state value given in Appendix B, one can get a lower 
limit for the time needed for a stationary flow to establish, 

(#o,obs. - #o,th.)/^ - \ohJ V * - 350 y- We therefore 
expect very little change in the morphology of the mid-IR or 
CO emission. The only appreciable evolution would involve 
the motion of the star relative to the IR or CO peak (1 " in 23 
yrs), since the shell is supposed to remain nearly static in the 
near future. 
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However, during its evolution towards a stationary dynam- 
ical regime, the velocity field must undergo a drastic varia- 
tion to reach the steady-state solution. CO mm observations, 
which provide excellent spectral resolution, might thus reveal 
significant velocity changes. A higher spatial resolution map of 
CO emission would be useful in this regard and of great help 
to better understand the kinematics underlying the remarkable 
pattern observed for line profiles in Fig. 1. In the next decade, 
ALMA will offer excellent opportunities for such observations. 
Numerical simulations of the time evolution of the wind/cloud 
interface (which, to our knowledge are not available) would 
also be very useful to indicate how the evolution of the velocity 
field will proceed in the early phase. 

In the above reasoning, we implicitly assumed that the 
dense shell will smoothly evolve towards the steady-state solu- 
tion but this is in no way evident. If indeed instabilities develop 
efficiently at the interface (as suggested by CH and CH + varia- 
tions), the cloud may simply be gradually destroyed as the star 
mov es. Su c h a picture would be consistent with the suggestion 
by iHerbigl (1 19581) that the absence of IS material south of the 
star is due to "clearing" along the path followed by HD 34078 
in the past (the clump involved in the second scenario might 
then simply be a fragment of the initial cloud in the process of 
photoevaporation) . 

Focusing now on the present state of HD 34078's envi- 
ronment, we note that it represents a remarkable PDR and 
shock for which many observational constraints are or might 
be available, thanks to the presence of a background UV- 
bright star. Geometrical parameters are now well determined 
and physical conditions in the ambient cloud relatively well 
constrained (sensitive CO emission observations further away 
from HD 34078 would allow us to better characterize them 
and verify that the peculiar excitation conditions determined 
in Sect. 2.2 are specific to the immediate vicinity of the star). 
Then, the HD 34078 PDR may be used as a reference to test 
our understanding of various physical and chemical processes 
occuring elsewhere at cloud interfaces subject to less extreme 
conditions. 
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Table 1. Observation parameters. The projection center of all the data is: 0:2000 = 05' l 16 m 18.15 s , 62000 — 34 deg 18'44.3". 
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Table 2. List of observations and measured CH and CH+ equivalent widths (in mA; uncertainties as estimated in Sect. 3.2.1 are 
given in upper index). a Resolution near CH A4300. 6 W43oo(CH) and W4232(CH+) values from OHP/ELODIE spectra have 
been corrected by +6% and +8% respectively (see text). c Column densities in units of 10 14 cm -2 , measured directly from line 
profiles for high resolution spectra (McDonald) or computed from W values using Eqs 1 and 2 for other data. 
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Appendix A: The contribution of the IC 405 
nebula to HD 34078 spectra 



B. 1. Shape and radius of the shell 



France et al. (2004) observed the diffuse emission from IC 405 
at four positions (their Posl to Pos4), two of which (Posl and 
Pos2) are located close to the HD 34078 position (75 " offset, E 
and W respectively). Thus, diffuse emission is certainly present 
towards HD 34078 itself and may contribute significantly to the 
flux collected in the FUSE apertures. To estimate this contribu- 
tion, we retrieved from the FUSE database the spectra obtained 
at Posl to Pos4 (cf. Fig. 1 from [France et al.l 12004). Pos2 is 



the brightest region; around 1050 A, the surface brightness is 
about 30% larger than at Posl, indicating spatial variations of 
this emission. Since Posl and Pos2 are symetrically located 
with respect to HD 34078, a first order estimate of the diffuse 
flux towards the star is the average of Posl and Pos2 spectra. 
We thus estimate that the diffuse flux received in the LWRS 
aperture is about 7% of the HD 34078 flux around A = 1050 
A (this fraction decreases with wavelength since the diffuse to 
ste llar flux ratio gets l ower at longer wavelengths, as discussed 
by lFrance et al]l2004l) . If the surface brightness is locally uni- 



form over the LWRS aperture centred on HD 34078 's position, 
the contribution of diffuse emission scales linearly with aper- 
ture size and should be about 1 1 times lower in MDRS spectra. 
Thus, while LWRS spectra of HD 34078 are significantly af- 
fected by diffuse emission, the MDRS spectrum should be es- 
sentially free of such "pollution", except possibly at the short- 
est wavelengths. 

To estimate the impact on HD 34078's LWRS spectra of the 
contamination by diffuse emission, one needs to examine the 
spectru m of the latter (the Posl spectrum is displayed in Fig. 
7 from lFrance et alJl2004h . Although its S/N ratio is limited, it 
is clear that it differs from that of HD 34078 in two respects. 
First, narrow lines appear to be fainter and shallower (only the 
strongest lines are detected); this is due, at least in part, to the 
lower effective resolution implied by the extented nature of the 
source. Second, broad H2 lines are narrower, indicating that the 
average pathlength of scattered photons through the molecular 
gas is smaller than that followed by direct HD 34078 photons. 
These two properties of diffuse emission qualitatively explain 
the peculiarities of the 8 th spectrum and indeed, one finds that 
by combining it and our estimate of the diffuse emission spec- 
trum towards HD 34078, it is possible to reproduce spectra 1 to 
7 fairly well. We conclude that the apparent changes in the 8 th 
spectrum can be attributed mainly to a smaller contribution of 
diffuse emission due to the use of MRDS instead of LWRS in 
the earlier spectra. 



Appendix B: The steady-state bow shock model 
confronted to observations 

Extensive work has been performed to describe the geom- 
etry and velocity structure of steady-state bow shocks (e.g. 
van Buren et aIlll990tlMa"c Low et al.|[l99lllWilkinll 1996b and 
in the thin-shell limit, there are simple analytical predictions 
that can be directly compared to observations. 



Assuming that the IR arc detected by lFrance et alJ (120071) cor- 
responds to a steady-state bow shock viewed in projection 
onto the sky, we can first check whether the apparent geo- 
metrical properties are consistent with the predicted ones. At 
first sight, the distance between the star and the apex of the 
bow shock looks too small as compared to its radius of curva- 
ture on the Spitzer image but we have to account for projec- 
tion effects which somewhat influence the appearence of the 
arc. HD 34078's tangential and radial components (Vj and Ir- 
respectively) are well constrained by observations; we adopt 
Vt = 100 km s _1 and V r = +59 km s _1 (in the LSR system; 
we remeasured the latter value from our own visible spectra). 
The velocity vector is then inclined by an angle 9 — 30° with 
respect to the plane of the sky (Fig.fPfl). 

To estimate the impact of projection effects on the h/R CtP 
ratio, wher e h is the apparen t standoff distance (h/D ~ 10 — 
20 " after iFrance etalJl2007l we adopt h/D = 15 " with 



D = 530 pc, the distance to HD 34078) and R c . p is the ra- 
dius of curvature of the arc as seen in projection on the sky 
(Rc,p/D ~ 37 "), we approximate the bow shock geometry 



by a paraboloid (Ivan Buren et al.l 1990: Wilkin 1996) and find 



after some algebra that 



, 3 + cos 2 9 n 
Acos 9 



and 



(B.l) 



R, 



2cos6» 



Ro, 



(B.2) 



where R is the distance between the star (S) and the apex 
(A) of the shock. With 9 = 30°, we get h = 1.08i? and 
h/R c ,p — 0.62 while the estimate inferred from the IR map 
is h/R c . p = 0.40 (the upper limit is 0.54 for h = 20 ")■ Thus, 
projection effects appear to be insufficient to explain the rel- 
ative large radius of curvature of the arc. Moreover, it seems 
difficult to explain in this mod el why the arc doe s not extend 
further southward (cf. Fig. 5 in lFrance et al.ll2007l) . Given that 
projection effects are in the end very limited, the standoff dis- 
tance can be estimated to be about 0.04 pc. 

Although the observed arc shape is not well fitted by the 
model prediction, the IR data can nevertheless be used to get 
a rough estimate of the distance between the star and the point 
where the line of sight intersects the shell (L2 in Fig. [I4l l. by 
making the reasonable assumption that the latter is axially sym- 
metric. From the 24 /im map, we estimate that d(S, L-^/D ~ 
40 ", implying d = 0.10 pc. Is this value compatible with the 
UV flux necessary to explain the amount of highly excited H2 
? B05 found that a radiation field about 10 4 larger than that in 
the local ISM is required, which is obtained at 0.2 pc from the 
HD 34078, a value in reasonable agreement with our estimate 
for d{S,L 2 ). 
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B.2. Momentum balance and ambient density 

In steady-state, the standoff distance, Rq, is set by a momen- 
tum balance equation (Eq. 2 in van Buren et al. 1990; Eq. 1 in 
Wilkin 1996) which is 



Rn — 



rh*V\ 



w 



(B.3) 



where m* is the mass loss rate, V s the terminal wind ve- 
locit y (10~ 95 Mq yr " 1 and 800 km s _1 respectively af- 
ter iMartins et a l. 2005), pa the mass density of the ambient 
medium and V s the star velocity (116 km s" 1 with the above 
values for Vj and Vr). The ambient H number density can be 
estimated either from C ( re = 700 cm" 3 ; B05) o r C 2 ab- 

sorption lines (re = 300 cm" . iFederman et al. I ll994l) . assum- 



ing the latter are not significantly contaminated by the dense 
shell. Adopting n — 500 cm" 3 , we get a "theoretical" value of 
i?o = 3.5 x 10" 4 pc while observations indicate Rq ~ 0.04 pc, 
in marked disagreement. In other words, an ambient density as 
low as of a few 10" 2 cm" 3 would be required for a steady-state 
bow shock to be at the observed distance, which is highly un- 
realistic for molecular-rich gas. Radiative pressure from stellar 
photons might help to maintain the shell at a distance larger 
than expected on the basis of the wind pressure alone. One can 
obtain easily an upper limit for the radiative to wind pressure 
ratio, P m r\/Pw, by assuming that all photons impinging on the 
shell are absorbed: P m( [/Pw = L/(rre*cVw). This ratio is of 
the order of 10" 3 for HD 34078; radiation pressure is therefore 
negligible here. 

Since P\y scales linearly with m*, one may wonder whether 
th e mass loss rat e has b een underestimated. Prior to the study 
bv lMartinseTa l. (2005), the adopted value for HD 34078 was 
10~ 6 ' 6 Mq yr _1 (i.e. larger by a factor of 800 than the present 
estimate !) and even with this much higher rate, the required 
ambient density would amount only to re ~ 20 cm" 3 . The 
much lower recent mass loss estimate is based on the availabil- 
ity of UV lines (C IV A 1550 mainly) which better probe weak 
winds; the uncertainty on the revised value is estimated to be of 
a factor of about 3 (F. Martins; private communication). Then, 
m* cannot have been underestimated by a factor large enough 
to explain the discrepancy between the observed and theoreti- 
cal steady-state R values. 

Stationary bow shock models also provide specific predic- 
tions for the mass surface density or equivalently the column 
density N(H) of swept-up material trapped in the bow shock (cf 
Eq. 7 from van Buren et al. 1990 or Eq. 12 from Wilkin 1996). 
iV(H) scales as rh 1 J 2 V^J 2 V^~ 1 n 1 / 2 and with the values quoted 
above, we get N(H) = 6.0 x 10 17 cm -2 . This prediction is to be 
compared to the H column density in the hot PDR component 
of B05: 2.7 x 10 19 cm" 2 , including H 2 only (their Tab. 4). Note 
that their exceedingly large predicted value for the H I column 
density was due to the assumption of steady-state equilibrium 
for the photodissociation of H 2 ; in our scenario, this assump- 
tion is no longer realistic. The observed value corresponds to 
the dense material along the line of sight to HD 34078 (i.e. lo- 
cated at L2 in Fig. [Pill while the "theoretical" one refers to the 
apex position (A in Fig. [14). This does not make a large dif- 
ference however since Wilkin's results (his Fig. 4) indicate that 



the surface density normal to the bow shock varies slowly with 
position away from the apex. One should also consider that in 
the geometry of Fig.[l4] the shell is not crossed normally by the 
line of sight but with an indication angle ip, but this involves a 
factor of at most a few. Obviously, this cannot explain the large 
discrepancy between the two values above. 

In the 24 /im map, the arc defines a roughly hemispherical 
cavity (with a radius of about R C , P ) and one can easily get an- 
other estimate for the column density in the compressed shell 
by simply assuming that material from the ambient cloud (i.e. 
with re = 500 cm" 3 ) initially filling this cavity has been swept 
by the stellar wind to form the present shell. This leads to a 
column density of N — rei? c p /3 and interestingly, this expres- 
sion provides a value, N ~ 5 x 10 19 cm -2 , comparable to the 
estimate of B05 for the hot PDR component. 

B.3. Velocity field of the compressed gas 

Another way to assess whether a steady-state bow-shock model 
is consistent with our data is to compare the observed and 
predicted velocity fields. Since 12 CO(2 — 1) emission traces 
dense gas within the shell of compressed gas, our CO data 
can be used to constrain the velocity field around HD 34078. 
One may wonder in particular, whether the double-peaked pro- 
files with their remarkable symmetry properties (Fig. 1) sim- 
ply arise from the fact that the IRAM beam intersects the 
paraboloidal wind/cloud interface twice. To compute a "syn- 
thetic" 12 C O(2 — 1) emissio n map, we adapted the model de- 
veloped by iPetv et al. I d2006l) to describe the outflow around 
HH30. We re lied on the analytical expressions provided by 
Wilkinl (119961) for the geometry, velocity field and mass sur- 
face density (see also Ivan Buren & Mac Low 1992 ), with the 
parameter values considered above. For practical reasons, the 
thickness of the boundary layer has been assumed to be Rq/20 
(while it is zero in Wilkin's model) and the absolute surface 
mass density has been scaled so as to reproduce the observed 
intensity. The underlying assumption is that the medium is op- 
tically thin which is reasonable given the strength of the CO 
emission. The resulting signal model was convolved with the 
IRAM-30m 230 GHz beam, an important step given the beam 
size relative to the source extent. Finally, in order to assess 
whether a specific model is acceptable or not, we compared the 
synthetized and observed maps of spectra as well as position- 
velocity and channel maps. 

Several difficulties arised when confronting the model with 
observations. First, the extent of the 12 CO(2 — 1) emission 
tends to be too limited if one adopts Rq/ D = 15 arcsec. This 
is related to the fact that the expected radius of curvature is too 
small for such a Rq value, as compared to the observed one 
(cf above). Since here we are mainly interested in the veloc- 
ity field, we simply adjusted Rq so as to match the observed 
extent of the emission. Second, the velocity range over which 
12 CO (2 — 1) emission appears in the model is much larger 
than the observed one. Indeed, the velocity field of the gas 
scales linearly with V s : in particular, the typical separation AV 
between the two CO emission peaks near HD 34078's position 
should be of the order of 0.7 V s while we observe only 0.02 V s 
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(corresponding to about 2 km s _1 ). Artificially modifying the 
star's velocity to the former value allows us to qualitatively re- 
produce the emission properties close to the star. However, the 
double-peaked character of the model profiles tend to be less 
pronounced that in the observed ones. Other velocity/density 
distributions might be considered to get a better fit, but clearly, 
only higher spatial resolution observations would allow us to 
get unambiguous constraints on such models. 

Another constraint can be obtained from optical absorption 
lines arising towards HD 34078. With the parameters quoted 
above, the expected velocity of the gas from the shell along the 
sightline to HD 34078 (i.e. at L 2 ) is Vlsr - -15 km s _1 . Such 
a shift between highly excited H 2 lines (tracing the dense layer) 
and absorption from species located beyond the shell would be 
easily detectable, but B05 failed to find any significant velocity 
difference between the two components. 

To summarize, the dense shell at the stellar wind/molecular 
cloud interface (with a density of about n ~ 10 4 cm~ 3 and 
column density N(H) ~ 3 x 10 19 cm~ 2 , thus corresponding 
to a thickness of ~ 10~ 3 pc) is located at a distance of the 
star that is consistent with the excitation of H2 but it does not 
display the properties expected for a steady-state bow shock: i) 
the arc is not curved enough; ii) the shell is too far from the star 
for the momentum balance to be satisfied; iii) the amount of 
material swept up by the wind is too large and iv) the velocity 
field shows very little deviation from the ambient value. The 
limited extent and somewhat irregular geometry of the arc are 
additional indications against a steady-state bow shock. 



